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ABSTRACT 

Biological net productivity, measured in terms of the change 
in biomass with time, affects global habitability and the quality 
of life through biochemical and hydrological cycles and by its 
effect on the overall energy balance (cf. Botkin (1982) or 
Botkin, et.al. (1984). Estimating leaf area for large ecosystems 
such as forests is one of the more important means of monitoring 
this productivity today. 

For a particular forest plot, the leaf area is often 
estimated by a two-stage process. In the first stage, known as 
"dimension analysis", a small number of trees are felled 
("sacrificed") so that their leaf areas can be measured as 
accurately as possible. The leaf areas of the sacrificed trees 
are then related to non-destructive, easily-measured features 
such as bole diameter or tree height, by using a regression 
model. In the second stage, the non-destructive features are 
measured for all or for a sample of trees in the plots and then 
used as input into the regression model to estimate the total 
leaf area. 

Because both stages of the estimation process are subject to 
error, it is difficult to evaluate the accuracy of the final plot 
leaf area estimates; indeed, it is not even possible to establish 
a meaningful criterion of accuracy without resolving issues 
involving sampling and modeling. This paper illustrates how a 
complete error analysis can be made, using an example from a 
study made on aspen trees in northern Minnesota. 
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1. BACKGROUND 

A joint study of leaf area was made in the Superior National 
Forest in northern Minnesota during the years 1983-4 by the 
National Aeronautics and Space Administration and the University 
of California at Santa Barbara. The main objective of this study 
known as COVER (Characterization of Vegetation with Remote 
Sensing) was to be able to relate the leaf area index (leaf area 
per unit ground area) of typical boreal .forest plots to 
remotely-sensed multispectral scanner data, with the ultimate 
goal of being able to estimate biological productivity by remote 
sensing. Botkin and Running (1984) provide a discussion of 
COVER'S basic approach. For a full description of the project 
the reader is referred to the forthcoming paper of Woods and 
Botkin (1986). 

For a particular forest plot, the leaf area is often 
estimated by a two-stage process. In the first stage, known as 
"dimension analysis", a small number of trees are felled 
("sacrificed") so that their leaf areas can be estimated as 
accurately as possible. In practice, leaf area for entire trees 
is not measured directly but is instead obtained from a ratio or 
regression model using leaf weight as an auxilliary variable, as 
in Watson (1937). The estimated leaf areas of the sacrificed 
trees are then related to non-destructive, easily-measured 
features such as bole diameter or tree height, by using a 
regression model. In the second stage, the non-destructive 
features are measured for all or for a sample of trees in the 
plots and then used as input into the regression model to 
estimate the total leaf area. 
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Because both stages of the estimation process involve 
errors, it is difficult to evaluate the accuracy of the final 
plot leaf area estimates. Some assessments of accuracy of 
dimension analysis regression fits have been found in the 
literature (cf. Pollard {1970 and 1972)), but they do not take 
into account errors in the leaf estimates of the sacrificed 
trees; in addition there do not appear to be any published 
estimates of accuracy at a plot level. In order to establish a 
meaningful criterion of accuracy, underlying sampling and 
modeling issues must first be resolved. For example, should 
finite-sampling or infinite-population models be assumed? How 
should model error variances be structured? Should sampling be 
random or should it be forced to cover the domain of the 
predictor variables? 

This paper describes how an extensive error analysis of the 
COVER leaf area index (LAI) estimates was made, with the intent 
of pointing out some of the practical problems that arose and how 
they were handled. Results given here deal only with plots of 
aspen trees, although stands of other species were also included 
in the study. 

The organization of the paper is as follows: Section 2 

describes in a procedural sense how the LAI estimates were made, 
both at the tree level and at the plot level. Section 3 provides 
the mathematical setting for obtaining the LAI estimates and also 
addresses the problem of accuracy estimation for both stages of 
the estimation process; Section 4 provides some numerical results 
including an error analysis, and the fifth and concluding section 
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discusses the strengths and weaknesses of the COVER procedure, 
suggesting how it might have been improved. 


2. THE COVER ESTIMATION PROCEDURE 
The target areas for which it was desired to estimate the 
LAI, consisted of 32 circular plots of (preponderantly) aspen 
trees. Each plot had a radius of 30 meters and was sampled by 
five circular subplots, (usually of radius 8 meters) arranged in 
a quincunx. The 30-meter plots were roughly the areas "seen" by 
a multispectral scanner carried aboard a helicopter hovering 
overhead (Figure 1). 

2.1 Leaf Area Estimation for Sacrificed Trees 

In COVER, the selection of the sacrificed trees was made so 
as to best represent the range of variation both in size and in 
local environment, of aspen trees in the region. Tree sizes were 
quantified in terms of "diameter at breast height" ( DBH ) , an 
ecological term which is defined to be 1/ir times the perimeter of 
the bole 1.4 meters from the ground. Sacrificed trees were 
chosen near but not in, the plots. They were selected at random, 
subject to two restrictions: (1) a minimum number of trees from 

each of a predefined set of DBH strata had to be filled; (2) 
unusually deformed or diseased trees were not taken, although no 
effort was made to restrict the sample to trees in the best of 
condition. 

The raw data for the sacrificed trees comprised various 
categories of measurements: those made for the entire tree. 
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those made for all branches, those made for selected branches, 
and those made only for small samples of foliage. Measurements 
made on the whole tree were non-destructive, consisting of DBH, 

H - the height of the tree and - the height above ground of' 
the point where the lowest live branch joins the bole 
(see Figure 2) . 

Measurements made on branches and foliage were destructive 
and were taken for each of three crown strata. Branches were 
partitioned into crown strata ("upper", "middle" and "lower") by 
first computing the depth of crown of the tree, defined as DOC = 
H-H^: then dividing the crown into three parts so that if h is 
the height at which a branch joins the bole, then the branch is 
in the lower stratum if h < H^+DOC/3; the middle stratum if 
+ DOC/3 <_ h < H-DOC/3; or the upper stratum if 
h ^ H - DOC/3. A leaf was considered to be in a particular 
stratum if its branch was in that stratum, regardless of the 
leaf's actual position in the crown. 

Given a tree whose crown was partitioned into strata as 
above, the folowing destructive raw data was taken: 

(1) Total leaf area for a sample of leaves from the i-th 
stratum (i = 1,2,3), where i = 1 refers to the upper 
stratum, i = 2 refers to the middle stratum and i = 3 
refers to the lower stratum. For a given leaf, "leaf 
area" is defined as its one-sided surface area. This 
was measured in the field by an optical scanner. 

(2) Total leaf fresh weight for each of the above samples. 
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(3) Branch dimensions - For each branch of the tree, (see 
Figure 3) measurements of: 

(a) D = diameter at the base of the branch 

(b) L = length of the branch (straight line distance 

from base to tip). 

(c) = straight line distance from the base of the 

branch to the point where the first live 
sub-branch joins the branch. 

(These measurements are analogous to DBH, H 
and H^, which pertain to the entire tree.) 

(4) Total foliage fresh weight for 

(a) a stratified random sample of branches for which a 
foliage-weight-vs-branch-dimension regression model is 
appropriate. These will be called "regression" 
branches in the remainder of this paper. 

(b) "special" branches for which a foliage-weight-vs- 
branch-dimens.ion regression model is not appropriate. 
Examples of these are "leaders’ - i.e. the very tops of 
trees, which contribute to the total leaf area, but are 
not true branches. 

We will refer to all branches in (4) above as "sampled" 
branches, an important subset of which is the collection of 
regression branches (4a). Only data from the regression branches 
was used to fit the model (see Section 3.1.1). 

A leaf area estimate for each sacrificed tree was made at 
the stratum level and then summed to obtain an estimate for the 
entire tree. For a given stratum the leaf area estimate was 
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obtained by applying an area-to-weight ratio to a leaf weight 
estimate using a regression model on branch dimensions. Because 
there were insufficient numbers of regression branches to allow 
separate models within strata, one model was used for all 
branches of each tree. The ratios were estimated by using a 
linear model to smooth sampling errors associated with raw ratios 
obtained from the data in (1) and (2) above. 

2.2 Plot-level Estimation - 

Once leaf area estimates were available on the sacrificed 
trees, they were regressed on functions of the non-destructive 
tree-level measurements to obtain a "dimension analysis" 
prediction model. During the second stage of the LAI estimation 
process, this model was then used on each plot to predict its 
total leaf area which was then normalized to a LAI estimate. 

3. MATHEMATICAL MODELS 

In order to construct the LAI estimates and a measure of 
their accuracies, some underlying conditions and response models 
must be assumed. This section discusses the problem of model 
selection and shows how the LAI estimators were obtained. We 
begin with the estimation of leaf area for the sacrificed trees, 
(within-tree models) then discuss the dimension analysis 
( between-tree ) model and finally, show how it is used to obtain 
the plot-level LAI estimates. While it would be possible to use 
one set of notation to identify relevant quantities in all stages 
of the estimation process, such notation would be 


6 


cumbersome. Consequently, the notation will be simplified, 
expanded or even changed according to the needs of each 
subsection below. 

3.1 Within Trees 


For a given sacrificed tree, let q^j be the total area of 
the foliage on the j-th branch in the i-th stratum; also let 
h = be the total leaf area of the tree and let q be its 

estimate. The main contributions to the error in estimating q 
come from the regression estimates of leaf weight and from 

a 

ratioing to leaf area. Whether q is biased or not and how its 
precision should be estimated depends on (a) the sampling 
strategy, (b) what models are assumed and (c) the form of the 
estimator. 

3.1.1 Foliage Weight Regression Model 

Let N be the total number of branches on the tree and let Y. 

3 

be the total leaf weight on the j-th branch, j=l,...,N (ignoring 
strata at present). The Y^ can be considered to have arisen in 
one of two ways: as a set of fixed numbers, or as realizations 

of a random process through time. Under the first formulation, a 
function of the Y's such as the total leaf weight, is estimated 
as a population parameter with the sampling design playing an 
important role. In the second case, the function is random and 
is thus predicted by a model-based approach. Further discussion 
of this topic may be found in Royall(1970, 1985), Kalton(1983) 
and Hansen, Madow and Tepping ( 1983 ) . 
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Preliminary studies indicated that a good auxiliary variable 
for predicting (estimating) Y^ from the available branch 
measurements D,L, and L^, was the branch depth of crown, 
dc = L-L^ . Plots of leaf weight vs dc for various trees 

2 

suggested the regression function E(Y^) = + f^dCj > 

where dc^ is the value of dc for the j-th branch. Since the 
foliage of a tree may vary throughout a year (or years) with 
little change in the branch dimensions it seems appropriate here 
to use the random ("superpopulation") model and think in terms of 
predicting, rather than estimating the Y's. 

Let Y = (Y 1 ,...,Y n ) T and let the N x 3 matrix 

T 2 T 

X = (Xi|x 2 l ** *I^N ) where = (If dCj, dc^ ) . A 

superpopulation model for leaf weight, conditional on the x^ was 

then initially taken as 

Y = X£ + £ (3.1) 

where ( 8 ^ 62 ' $ 3 ) and £. i s an error vector. 

It was observed (see Figures 4a-c, for example) that (3.1) 
fit more accurately for small branches than for large. 

Examination of the pooled residuals from preliminary regressions 
of Y on X for all 32 trees showed that their variance was roughly 
proportional to the square of dc; hence it was decided to 


represent the variance of £ by 

v ar(_c|X) = a 2 V, (3.2) 

2 

where o is an unknown scalar, and V = diag ( v i • v o' * * * v n) 



Empirical models such as (3.1) are only intended to 
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represent a local relationship; hence in order for one to have 
confidence in the regression predictor of Y, not only should the 
fit of the model to the regression data be good, but also 
prediction should be avoided for x-values outside the range of 
the regression data. 

Figures 4a-c which plot leaf weight against dc for trees 
#21, #24 and #27 respectively, illustrate typical situations 
which can occur. In these figures, vertical bars have been added 
above the dc-axis showing the distribution of dc values for the 
unsampled branches. In Figure 4c (Tree 27) the fit is good and 
the dc values for all the unsampled branches lie well within the 
domain of the regression points. In Figure 4a (Tree 21), the fit 
is not very good, but reasonable predictions can still be made 
for the unsampled branches because their dc-values lie within 
those of the regression branches. In Figure 4b (Tree 24), the 
fit is not bad, but there are three large branches with dc-values 
outside the domain of the regression, thus requiring 
extrapolation. In a case such as this it would have been better 
to not sample the branches randomly, but instead to have made 
sure that the largest branches were included in the regression, 
as suggested by Royall (1970). 

Let I_* = diag (c 1 ,..,c N ) where c^ = 1 if the j-th branch is 
a regression branch, otherwise c^ = 0. Standard estimation of 3^ 
with the model given by (3.1) and (3.2) consists of finding the 

vector b which minimizes L(b) = ( Y-X b) T Y. * Although 

the b obtained in this manner is unbiased and has desirable 
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necessarily produce useful predictions of Y for all dc values 

likely to be encountered. In particular, foliage weight 
T 

predictions b x . were sometimes negative for small branches, and 

in one case, where a tree had many small branches, the estimated 

foliage weight for the entire tree was negative. To prevent this 

anomaly, the model was revised to constrain £ to be nonnegative 

and b was then calculated as the solution to minimizing L(b) 

subject to b > 0. It is known - e.g. see Lawson and Hanson 

(1974), that the constrained b is the solution to one of the 

unconstrained subset regressions using some or all of the inde- 

0 2 

pendent variables dc -1, dc and dc ; in fact it is the solution 
to the subset regression which has the smallest value of L, given 
that the (unconstrained) regression coefficients are nonnegative. 

Even when e is assumed to have a multivariate normal 
distribution, the distribution of b, or even its variance, is not 
known. In order to obtain an approximate expression for Var(b), 
we pretend that the model (3.1) originally included only those 
independent variables whose subset regression produced the 
optimal nonnegative b, so that x ^ , J3 and b are now understood to 
be p-dimensional vectors and X is a N x p matrix, where 
1 < p < 3, depending on the tree. In terms of the revised X, we 
have 

b = (xVx)'Wy (3.3) 

where W = 1^ V * . 

3.1.2 The Ratio Model 

In addition to the above model for describing the 
relationship between leaf weight and branch dimensions, another 
is required to relate leaf area to weight. 
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as aescrioed in Section 2.1, samples of foliage were taken 


from each stratum of each tree to obtain leaf area and leaf 
weight measurements. Let y^ k an< ^ x ^ denote the respective total 
area and total fresh weight of the leaves in the sample from the 
i-th stratum in the k-th tree and let and X^ k be the 

corresponding population quantities (i=l,2,3; k=l,..,K). 

To a first approximation, it can be assumed that the sample 
ratio r^ = y^/x^ is an unbiased estimate of the population 


ratio R ik = li lk /X ik ! i.e. 

r ik ' R ik + e ik 


(3.4) 


where e ik is an error term. In general, one would expect r^ k to 
be a more accurate estimate of R^ k when the sample is large? i.e. 

when x^ k is large. For this reason we assume Var (^^Ix^) = 

2 2 
• (0 /x^ k for some positive constant w . 

Rather than base the estimation of leaf area of the 
sacrificed trees on the rik directly, it was decided to first 
smooth the raw ratios using an additive model for the true ratios 
which allow one to. incorporate information about all sacrificed 
trees and strata to obtain the estimate of the ratio for a given 
tree and stratum. The true-ratio model is given by 


R ik = v+ “i + 6 k + d ik 


(3.5) 


where y is an overall mean, is a stratum effect, 3 k is a tree 

effect withZa^ = = 0, and the d^ k are mutually independent 

2 

random errors with E(d^ k ) = 0 and Var(d^ k ) = x • 

Let p^ k = E(R^ k ) = y+ a ^ + 3^. The additive model for the 
p^ k is partially justifiable by ecological considerations. It is 
a recognized fact (cf. Talbert and Holch (1957)) that for a given 



tree, "sun" leaves (those high in the tree crown) tend to be 
thicker and smaller in area than "shade" leaves (those low in the 
crown) ; thus creating a tendency for the area-to-weight ratios to 
be ordered by strata with the lowest ratio for the upper stratum, 
and the highest ratio for the lower stratum. If a ^ < a 2 < 013 , 
equation (3.5) imposes this theoretical ordering on the expected 
ratios for each tree. 

Combining (3.4) and (3.5) gives 

r ik =y+a i + 3 Ic + e' i u (3.6) 


ik 


where e'^ k = + d^ k . It * s not unreasonable to assume the e^ k 


and d^ k to be uncorrelated in which case 

v a r <e' lk |x i( .) - w 2 /x ik + t 2 


(3.7) 


Let p^ k = y + + 3 k where y, and 3 k are obtained 

from fitting the model (3.6) by weighted least squares with 
weights proportional to x^ k . It is advantageous (in terms of a 
smaller prediction error) to predict R ik by r ik when the second 
term of (3.7) dominates, or by p^ k when the first term dominates. 
In this study, we assumed the latter case to be true because (a) 
the sample sizes of leaves were fairly small, and (b) it was felt 
' that the R . , were well represented by the additive model - i ,e . 

X K 

2 

T was relatively small. 

Let p be a 3K x 1 vector of the p^ k and let r* be a 
corresponding vector whose components are the r^. Then one may 

** it ic 

write p = £ _r , where U is a function of the design matrix of 
the model (3.6) and the x.^. For purposes of computing the 
estimated mean squared prediction error (MSPE) of ri for a given 
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tree (say the k-th), it is assumed that the 3x1 vector p k = 
T 


(P 


lk 


, P 2 ](' p 3k^ i- s d i str i buted as multivariate normal with mean 


— k = ^ P lk' p 2k' p 3k' and covariance matrix y k , where y k is the 
submatrix corresponding to the k-th tree, taken from the 


covariance matrix y of p . The matrix _y_ is of the form 
2 2 

w A^ + x A 2' wbere and A 2 are f unct i° ns of £ and the x^; 

★ 

however V cannot be estimated unbiasedly because there is no 

2 2 2 2 
way of estimating w + t A 2 directly or m and x separately 

with the available data. An approximately unbiased estimate of 

y* is y* = £ 2 U D _1 (x)U T , \ where D(x) = diag(x^ k ) and w 2 is the 

residual mean square given by 

» 2 = UTx ik (r ik - P ik ) 2 J /2 (k-1 ) (3.8) 

2 

This estimate would be unbiased if x were zero. 


3.1.3 The Estimator of n and its MSPE 

Reverting to the discussion of the prediction of n for a 
given tree, we drop the subscript "k" of Section 3.1.2, using p 
to mean the 3x1 vector of smoothed ratio estimates and _y_ to 

a 

mean the 3x3 covariance matrix of p for the tree. Likewise, 

"R" will designate the 3x1 vector of true ratios for the tree 
with E (R) = p. 

The foliage weight estimate for a stratum was the sum of two 
components: one being the sum of measured weights for sampled 

branches; the other being the sum of predicted weights for the 
unsampled branches using the model (3.1). Let e. be the N x 1 
membership vector for the sampled branches of the i-th stratum; 
i.e. the v -th element of e^ s is zero or one, being unity if and 
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a sampled branch. Similarly, let be the membership vector* 

for the unsampled branches in the i-th stratum. Then the 
estimated foliage weight for the stratum is 

^ m rri 

y . = e Y + e. X b 
J 1 —is — iu— — 

= (e T + ef P W*)Y (3.9) 

—is — iu — 

where P = X(X T W X) _ 1 X T . Let Z = (Yj/ V 2' Y 3 )T * The total leaf 

area estimate is given by 

A Afp A 

n = z p (3.10) 

In terms of y^, the actual foliage weight for the i-th 

T T T 

stratum can be expressed by = (e + e ^ u )Y = Y' where 

T 

e. = e. + e. . Let = (y,, y ? , y-,) . The actual leaf area 

“1 “1 S ™*X U X Z j 

T ~ 

is thus equal to £ 5.* To compute the MSPE of h , let the 12x1 
vector z be defined by 


* = <?, i T . P T , r t ) t 


and the 12 x 12 matrix A be given by 


(3.11) 



(3.12) 


ip ^ ^ 

where 1^ is a 3 x 3 identity matrix. Then z A z = y_ p_- R; i.e. 

a 

it is the prediction error n ~ n . 

Let V_ = E (y ) = (e* X _B_, e J X J, £3 X &) T - From the revised 
version of (3.1) and from (3.9), it can be seen that 
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E(^) = E(^) = ]j_. From previous assumptions we also have 
E (p ) = E(R) = £. Assuming further that y_ is independent of p 

rr» rp rr» ^ 

yields E(z A z) = p__P - p _P = 0, so that the MSPE of n is the 
T 

variance of z A z. 


Let E(^) = £ and Var(z^) = Z. If £ is assumed to be 

distributed multivariate normal, then 

MSPE(n) = Var(z T A z)= 4^A S A ? + 2 tr(Z A E A) \ (3.13) 

^ ^ 

We have already seen that ? - (J£ rP. '£ • £. ) • To evaluate 
(3.13) it remains to specify E which is of the form 



Z = 



ll2 

h.2 

0 

0 


0 



0 



(3.14) 


where E^, £ 22 and —12 are x 3 matrices Var (y_) , Var (£) and 




iji 

Cov (£,£) , respectively, and ~ — 12 * In terms of (3.12) and 

(3.14), equation (3.13) becomes 

MSPE(n) - u T Pf + T 2 (i.-u-u T )] lL + a T (J 11 -z 12 -j 21 + l22 )p 
+ tr[ Lll i + t 2 ( S . 22 -0 E 21 - U T I 12 )I (3.15) 

The elements of 1^,, Z 12 and l 2 2 are derived from the model 

a 

given by (3.1) and (3.2), and from the definition of in (3.9). 


15 


Specifically, it can be shown that 


Cov(y.,y.'') = a 2 [e^ s V e.,,. + e^P e . , r +e^P e^+e^P e-.J (3.16) 


* * 

where e. = I e- = I_ e. , It can be similarly shown that 
— i r ■“ is “"1 

Cov (y f -, y v ) = o 2 [£i s Z 2l , s + e* u P e., r l 


(3.17) 


Cov (y i ,y i i- )= a 2 eTv e £ , 


(3.18) 


The MSPE of n is estimated by substituting estimates for the 
unknown parameters in (3.15); i.e. ^ for y, p for p_, Y for and 

a 

. for I_. . (i,j = 1,2). The matrix V is the appropriate 

^ if ** 

submatrix of of Section 3.1.2, while the are estimated by 

"2 2 ~2 
substituting a for a in (3.12a-c), where a is the residual 


mean square given by 

a 2 = (Y - X b) T W*(Y - X b)/(n-p) , 


(3.19) 


and n = is the number of observations in the foliage weight 

2 

regression. In this analysis, terms of (3.15) involving x are 

2 

neglected; however* (1) x is presumed small and (2) the effect 
2 2 

of T is absorbed into £5 and hence y . Ignoring the other blocks 

2 

of £ involving x is equivalent to treating R as a fixed, rather 

than a random quantity, thus making (3.15) more like a "variance" 

rather than a MSPE as far as the ratios are concerned. If 

multiple observations of area-to-weight ratios had been taken 

2 

within strata, it would have been possible to estimate x . It is 
suggested that this be done in future studies of this type. 

3.2 Between Trees - The Dimension Analysis Model 

Given estimates of leaf area for the sacrificed trees, the 
next step in the LAI estimation, process is to relate these 
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estimates to the nondestructive measurements L, L^ and DBH. 

Let x be a p x 1 vector of functions of L, , and DBH which 

predict leaf area, n through a model 

n = Y T x + d (3.20) 

where d is an error term, independent of x, with a zero mean and 

a variance v(x) (conditional on x) which is an increasing 

T 

function of E(n ) = Y x. (The actual functions of L, and DBH 
producing the components of x are described in Section 4.2.) The 
leaf area for a new tree with given x-vector x' / may then be 
predicted by 

n = £ T x' (3.21) 

where £ is an estimate of Y obtained from the sacrificed tree 
data. 

3.2.1 Estimation of Model Coefficients 

Let n^, x^ and d^ play the respective roles of n , x and d 

in the model (3.20) for the k-th sacrificed tree (k = 1,..,K). 

Also for the k-th sacrificed tree, let n k be given by (3.10). To 

* ^ 

a first approximation is an unbiased predictor of ^ with MSPE 
given by (3.15). If the within-tree estimation errors are 
independent of the between-tree errors d^,/ then the 
^ follow a revised model 

\=I T S k + < (3.22) 

» 

with d^' being similar to d^, except that Var(d k |x. k ) = vfx^) 

+ Vq^, where v^ k is the within-tree MSPE given by (3.15) for the 
k-th tree. 

Ideally, it would be best to use the weighted least-squares 
estimate g = (X T V~ 1 X ) -1 X T V _1 n for g in (3.21), 
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= i > n 2 / • • »n R ; / a = 1 2£o I * * * l-K ) and 

V = diag £v(x k ) + Vq^]; but the function v(.) is not known and 

can only be roughly estimated from the data. Unlike the 

within-tree situation for which there were many observations of 

residuals available for estimating a weighting function, there 

were only 32 here. With such a small sample size, estimating the 

weights from the data could seriously bias the estimate of Y. As 

a consequence, the unweighted quantity 

T -1 T'* 

2 = (XX) X n (3.23) 

was used as an estimate of Y. Of course, cj is still an unbiased 

estimator of Y under the model (3.22). Its covariance matrix is 

given by T _ 1 T t l 

W = (XX) XV X (X X) . (3.24) 

3.2.2 Estimation of the Model Variance Function 

For purposes of estimating the precision of the plot LAI 

estimate, it was necessary to have at least rough estimates of 

the function v(.) and the matrix W. These estimates were 

obtained through a parametric model using the residuals of the 

unweighted regression as observations. Let e_ = fi - X £ be the 

K x 1 vector of residuals. Under the model (3.22) the covariance 

matrix of £ is 

Var (e |X ) = V - X (X T X) -1 X T V - V X (X T X)“ 1 X T 

+ X(X T X) _1 (X T V X) (X T X)“ 1 X T • (3.25) 

T 

and the residual sum of squares e e_ has expectation 

E(e T e) = tr [V] - tr [ ( X T X ) -1 (X T V X ) } . (3.26) 

The non-leading terms of both (3.25) and (3.26) are of order 1/K. 
Neglecting them is equivalent to treating the residuals as actual 
observations of the d k ’ in (3.22). Assuming that this is a 
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(3.27) 


reasonable approximation gives 

T* 

E(e ®)“ Ev(x k ) + Zv 
k k 

and 

Var(e R )= v ( x k ) + v Qk (3.28) 

where e k is the k-th component of e . 

Suppose now that v(x) can be modeled by 

v(x) = 9 0 (Y T x) 9 1' (3.29) 

where 0 q and 0 ^ are parameters. From (3.27), we have E(e^e)= 

. y q - * 

e 0 ^ 1 + j^Ok* G1ven an estimate of ®l» say 0 X , 0 O can be estimated 
by setting e T e equal to its estimated expectation; i.e. 

®n = (£ T £ - 2v 0 k)/ 2 (£ T x k ) ®1 (3.30) 

k k 

= GCBp , say, 

where 9g is the estimate of 0 g. 

Consider the normalized residuals 
z k = e k /[G(0 1 ) (£ T x k )°l ] 1/ ^ 2 , (k=l,..,K) where 0 ^^ is a trial value 
of 0^. Given the assumptions of this section, if 0 ^ is close to 
0 ^, the z k should have approximately unit variance; in fact, if 
the residuals are normally distributed, the z k should resemble a 
random sample from a N (0 , 1 ) population. We therefore choose 0 ^ 
to be the value of 0 ^ which minimizes the Kolmogorov-Smirnov 
statistic D, which is the maximum absolute deviation between the 
empirical distribution function of the z k and that of a standard 
normal distribution. Once 0^ is found, q q is already available 
as G( 0 g). Let v Qk be the sample analog of Vg k » For an arbitrary 
x, v(x) is then estimated by 
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0.31) 


V (x) = e 0 (2 r 2i) 01 

and W is estimated by using v = diag ( v(x k ) +v ok * ' fo2r Y. ^ 

(3.24). 

3.3 The LAI Estimate and its Precision 

Given £, the total leaf area of all trees in thm- sudoplotts 
within a plot can be estimated by the application of: 0,21)} since 
x is known for these trees. The total leaf area is: thenv 
normalized by the total subplot area to obtain the iL'J&fc esttiima fce . 
There are two main independent contributions of error lira: this 
estimate. The first is a between-subplot sampling tr.rosr dtae to 
the nonhomogeneity of the LAI throughout the plot. The: aecxarad 
error is the discrepancy between the model-predicted* I*,eax5: aunea 
and the actual leaf area within the subplots. It is not 
unreasonable to assume these two errors independent, iru which 
case the overall MSPE of the LAI estimate is 

E( A - A ) 2 = E(A-A) 2 + E ( X -X ) 2 (13. 32) 

where A is the estimated LAI over the subplots (and; aJs^tc is the 
LAI estimate for the whole plot),. A is the true LAI cs£ this plot 
and A is the true LAI over the subplots. 

For purposes of evaluating the first term of ({3)^02.)) „ the 
leaf area pattern for the whole plot is considered liLasedi s<n> that 
the LAI’s of the subplots become fixed numbers once fch.es siiisibplots 
are chosen. For the second term, we shall conditic.ni on- the 

observed x-values for both the sacrificed and the f.A<ofc trees. We 
thus use 
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Q = EH X - A ) 2 jLI + El(\ - A ) 2 | x j 

= Q x +0 2 (3.33) 

instead of (3.32) as a measure of precision of X. Here, 

denotes conditioning on a fixed plot leaf area pattern. 
3.3.1 Between-Subplot Sampling Variance 

The variance of the between-subplot sampling error is 
difficult to estimate without knowledge of a leaf area intensity 
function over the whole plot? i.e. the leaf area "pattern”. As a 
practical alternative, we pretend that the whole plot consists of 
M subplots, of which we have selected m at random. Let f denote 
the sampling fraction m/M and let and denote the respective 
true LAI and area of the i-th subplot (i = 1,..M). For the i-th 
sampled subplot let and a^ be the respective LAI and area 
(i = l,..,m), and let X^ be the estimate of \ . obtained from the 
dimension analysis model. 

We shall now condition on the leaf area pattern for the plot 
so that the A^ may be regarded as fixed numbers, as opposed to 
the random variables they would be under the model (3.20). With 
this restriction, the X^ are random only because of the subplot 
sampling process. The LAI over the m sampled subplots is then 

m m 

A = y(a./a*)A. , where a. = la. . For the whole plot, the LAI is 

A = j(A^/A.)A ^ where A. = % A^ . The difference X-A is the error 

due to sampling the plot. When the areas of all M subplots are 

equal, X is an unbiased estimator of A; however when the areas are 

unequal, X becomes a ratio estimator and is only approximately 

2 

unbiased. Assuming the bias is neglibible, = E{( x ~ A) |LJ = 
Var [X|LJ. In the general case of unequal A^, we have the 
approximation (cf. Cochran, 1977, p.155) 
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/ IVI , 

q= [(l-f)/ra]M^ E (A /A.) /(A. "A; x , 

i i 1 1 

If the A^ were known, Q-^ would be estimated by 

q x = (1 - f)m j - X) 2 /(m-l) (3.35) 

where f = a. /A. replaces f in (3.34). Since the A. are not known, 

we replace them with their estimates A^ to obtain 

Q = (l-f)m ?(a./a. ) 2 ( X ± ~X ) 2 /(m-l ) (3.36) 

1 i 

as the estimate of Q^. 

3.3.2 Within-Subplot Prediction Error 

Let be the true leaf area of the j-th tree in the i-th 

subplot and let ^ be the vector of leaf area predictors for 

that tree, with n** = ZZn ± ^ and x. . = ELx^. The estimate of the 

~ T 

LAI for the m sampled subplots is A = £ x../a. and the actual LAI 
is A = ^(a^a.) A^ = n../a. From (3.20), n^j = + d ij where 


E[dij|Xij] = 0, and Varfd^Jx^l = vtx^). Since £ is an 
unbiased estimator of y, the conditional mean prediction error of 

a 

A - A is zero; i.e. 


E(A -A|x) = E(£ T x. ./a. |x)-E(h . ./a. |x) (3.37) 

= (y T x. .— Y T x. . )/a. = 0 

The second term (Q 2 ) in (3.33) is thus the conditonal variance of 
A - so that 

Q 2 = E f ( A -A ) 2 ! x I = Varl (A -A ) |x] 

T 

= Var$(£ x../a. - |j../a.)|xJ 


= (1/a. 2 ) [xT.Wx. .+ 

Its estimate, Q 2 is obtained by direct substitution; 

Gj = ( 1/a? ) [xT.W x ..+ ? s .v( x ..)] 

z — — — . l j i j 

where W and v(.) are the respective estimates of the 


(3.38) 

i.e. 

(3.39) 

matrix W and 
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the function v(.) from Section 3.2.2.. The term x. .W x . . reflects 
the fact that Y.' is estimated, while the other term ZZv(x. j) 
reflects the actual model error. By increasing the number of 
sacrificed trees, the former term can be made arbitrarily small, 
but not the latter. In Section 4.3, the effect of each of these 
terms on Q 2 is evaluated numerically to give a feeling for 
whether the sample size of sacrificed trees was adequate. 

4.0 NUMERICAL RESULTS 

In this section, numerical results from the COVER experiment 
are given. For the most part these results are organized along 
the lines of the theoretical development in Section 3 so that a 
subsection numbered "4.x...y'' consists of a numerical 
illustration of the topic covered under Section 3.x...y. 

4.1 Within Trees 

Recall that in the COVER experiment 32 trees were 
destructively sampled so as to best estimate their leaf area; 
that for a given tree, foliage weights were estimated with the 
aid of a regression model and then ratioed to leaf areas for each 
of three crown strata. These trees are identified by a number 
ranging from 1 to 32 which corresponds to the order of their 
selection. Table 1 lists the trees in order of size (i.e. DBH ) 
showing the Tree Number, the DBH in cm, H in meters, the total 
number of branches (N), the number of sampled branches (NS) and 
the foliage weight regression sample size (n). 

4.1.1 Foliage Weight Regressions 

From Table 1 it can be seen that the foliage weight was 
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measured for anywhere between 20 and 100 percent of the branches.' 
For purposes of measuring foliage weights, it was originally 
intended that a census of branches on the sacrificed trees be 
taken. For the first five trees, a census or near-census was 
made, but the excessive amount of work encountered soon made it 
clear that some sort of sampling was necessary. For this reason 
the sampling fractions for the other trees were considerably 
smaller (except Tree 11 which had only 9 branches, all of which 
were sampled). 

While the relationship between leaf weight and branch 

dimensions was not overly strong, it was usually good enough to 

justify the use of regression to estimate the foliage weight 

totals within strata for trees in which branches were sampled. 

For each tree, Table 2 gives the dimension (p) of the revised 

regression model of Section 3.1, the regression coefficients 

given by (3.3) and the t-values associated with the coefficients, 

assuming the p-dimensional model. The coefficients b^ are scaled 

to obtain the predicted foliage weight (y) in grams, from the 

branch depth-of-crown (dc) in meters, from the equation 
** 2 

y = b Q + bj(dc)+ b 2 (dc) . Terms deleted from the model by the 
constrained least-squares algorithm are indicated by a dash " — 

In only 7 of the 32 regressions were all three possible terms 
included in the revised model, while two terms were used in 24 
cases and one term in one case. 

For many trees, more terms could have been deleted because 
they did not contribute significantly to the regression fit; 
however for simplicity, any terms with nonzero coefficients as 
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4 . 


found by the algorithm, were retained. From Table 2 it can be 

seen that in most cases where the intercept term was included in 

a model, it was found to have a small t-value, but that the 

linear term was usually significant and the quadratic term was 

almost always significant when included in a revised model. 

4.1.2 Estimation of Area-to-Weight Ratios 

The average raw ratio of leaf area to leaf weight was about 
2 

4.9 m /kg. As predicted, these ratios tended to be greatest in 

the lower crown stratum and smallest in the upper crown stratum, 

although the magnitude of the difference was fairly small 

(generally about 10 percent). For 26 of the 32 trees, the 

upper-stratum ratio was the lowest of the three, while in 21 

cases the lower stratum ratio was the highest. For 19 trees, the 

three ratios ranked exactly as expected; i.e. inversely with the 

crown position. The tree effect was even more noticeable, with 

2 

the average raw ratio ranging from about 4.3 m /kg for Tree 29 
2 

to 7.3 m /kg for Tree 30. In general, the largest ratios were 
observed on small trees. 

As explained in Section 3.1.2, it was decided to replace the 
raw ratios (r^ k ) by smoothed ones (P^ k ) which were felt to be 
better predictors of the actual ratios. Figure 5 shows the r^ k 

a 

plotted against the p^ k for i = 1, 2, 3; k = 1,2,..., 32. Points 
lying between the two lines are those cases (73 of 96) for which 


p.. was within 5% of r .. . Without replicated samples, it was not 


possible to formally test whether the additive model (3.6) was 


correct, but it appeared to fit the data well. The ot ^ were in 
the "correct" biological order (decreasing with i), thus so were 
the for each k. Estimated coefficients of variation of the 
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ranged from about 4 to 12 percent. It was not possible to 
calculate a comparable quantity for the raw ratios. 

4.1.3 Leaf Area Estimates for the Sacrificed Trees 

For each sacrificed tree, (3.10) was used to obtain n the 
estimated leaf area, and the sample analog of (3.15) was used to 
estimate its MSPE. For the sacrificed trees arranged in order of 

a 

increasing DBH, Table 3 gives r) , its standard error SE = 

(MSPE)^ 2 and proportional contributions to the MSPE due to 
branch sampling and to area-to-weight ratio estimation. These 
components were obtained by evaluating (3.15) with certain 
changes to _E* Specifically, if a census of branches were taken, 
thus making the leaf weights of all branches known (assuming 
negligible measurement errors), then in the expresion for £ given 
^ and would both be equal to Z_22 because y would 

be equal to y. As a result, the middle term of (3.15) would 
vanish. Let MSPE1 denote (3.15) evaluated with!, changed in this 
way. The normalized difference (MSPE-MSPEl )/MSPE is then shown 
in Table 3 under the heading "BRANCH SAMPLING" as the 
proportional contribution of branch sampling to the MSPE. 

Similarly, if the ratios of leaf area to leaf weight were 
already known on the stratum level, one would replace £ with R 
so that in Z, would be replaced by x 2 X anc * U=I. Let MSPE2 

i 

denote (3.15) evaluated with V = x 2 I and U=I. The first term vanishes 
regardless of the value of t 2 . We have also neglected the third term under the 
assumption that x 2 is negligible . The proportion of MSPE due to ratio estimation 


by (3.14) , ^ 
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then equal to 100(MSPE - MSPE2)/MSPE and is shown in Table 3 in 
the last column entitled "RATIO ESTIMATION”. A perusal of the 
two reductions (ignoring those trees for which a census had 
already been taken) , reveals that neither was consistently 
greater , thus indicating that the manpower expended was 
well-allocated between the tasks of collecting total foliage 
weights for the sampled branches and measuring areas and weights 
for the small foliage samples. 

To quantify the effect on the overall plot LAI estimates, of 
errors in the leaf areas of the sacrificed trees, the MSPE of the 
plot estimates was evaluated with and without the effect of these 
errors, and then compared. (See Section 4.3). 

4.2 Between Trees 

4.2.1 Fitting the Dimension Analysis Model 

Previous work on non-destructive dimension analysis models 

for aspen trees (cf. Peterson, et al (1970) or Pollard (1970 and 

% 

1972)) has utilized DBH as the main predictor of leaf area. The 

additional inclusion of tree height (H) to linear or loglinear 

models was at best marginally better. For this study. Woods and 

Botkin (1986) obtained the measurement which made it possible 

to use information about the size of the tree crown, say its 

volume, as a predictor. This was done by modeling the crown as a 

cylinder or cone which enables one to approximate the crown 

2 

volume by c(H-H^)w , where w is the crown width at its base and c 
is a constant. The crown width was not directly measured in this 
study, so DBH was used as a surrogate. As a result, functions of 
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the dimension analysis model. Figure 6 illustrates the 

relationship between leaf area and VOC by plotting the former 
1 /2 

against (VOC) for the 32 sacrificed trees. 

The DBH values for the trees encountered in the subplots 

'ranged from less than 1 cm to about 40 cm. Preliminary studies 

1/2 

indicated that the leaf area was roughly proportional to (VOC) 

2 

for small trees and to (VOC) for large ones. Rather than suffer 

the awkwardness of working with two separate models, it was 

decided to use both powers of VOC as predictors, so that in 

1/2 2 T 

(3.20) the vector x = [VOC , VOC ] . The prediction equation 

(3.21) thus is 

h * g x (VOC) 1/2 + g 2 (VOC) 2 (4.1) 

a 

where n is the predicted leaf area and the coefficients g^ and g 2 

2 

are given by (3.23). For leaf area measured in m , H and m 
m, and DBH in cm, g-^ and g 2 were calculated as 0.39476 and 
2.94887 x 10“ 7 respectively. 

4.2.2 The Error Variance Function 

Using the methods of Section 3.2.2, the estimated variance 
function v(x) given by (3.31) was calculated as 

v(x) = .35858 (n) 1 * 4 (4.2) 

Recall that v(x) estimates v(x), the conditional variance of a 
tree's leaf area given x. Taking values of x. observed on the 
sacrificed trees as typical, the estimated coefficient of 
variation [v(x)] ' /r\ ranges from about 0.6 for a small tree (DBH 
about 2 cm.) to about 0.15 for the largest trees (DBH about 
35 cm . ) . 
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4.3 Plot-Level Results - Error Analysis 

The leaf area index was estimated for the 32 plots using 
(4.1) to obtain an estimated leaf area for each tree in the 
subplots, then summing and normalizing by the sum of the subplot 
areas. Standard errors were estimated as [Q 1 + Q 2 ] 1//2 where Q 1 
and Q 2 are given by (3.36) and (3.39) respectively. Table 4 
shows the sampling fraction (f), estimated LAI (LAI), its 
standard error (SE) and its coefficient of variation (CV) for 
each plot. In this table, the plots have been ordered by their 
LAI estimates. 

^ «•* 

By looking at Q 1 and 0 2 separately it can be seen whether 
large standard errors in LAI can be attributed to the 
heterogeneity of a given plot (as quantified by ) or to the use 
of the dimension analysis model (as quantified by 0 2 )* A further 
breakdown can be made by separating the two terms of Q 2 * The 

~ / 2 ) —2 T ~ 

first term, Q 2 = (a.) x**H2£**' reflects the uncertainty in 


the estimates of the regression coefficients and could be made 

arbitrarily small by taking a sufficiently large sample of 

sacrificed trees. The second term = (a.) £Zv(x. .), 

•<£ l j 

represents the true model error and/or lack of fit. This term 


cannot be reduced by increasing the number of sacrificed trees 


In Table 4, the three components of the estimated MSPE - 
~ ( 1 ) ~ ( 2 ) 

02 ' °2 and 0^ are shown for each plot under the respective 

headings "COEF EST VAR", "MODEL VAR" and "SAMP VAR". In every 
case, the dimension analysis model was by far the smallest 
contributor to the MSPE. The ratio of the sampling component 
(Q^) to the coefficient estimation component (Q 2 ^), is shown 
under the column headed "RATIO". An examination of this ratio 
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ror zu or tne 32 plots, the coefficient estimation * 
component exceeded the sampling component, thus suggesting that < 
it would have been worthwhile to estimate the regression 
coefficients more accurately by taking a larger sample of 
sacrificed trees. 

In a few plots (e.g. Nos. 20, 21 and 89) the sampling 
component especially predominated. An examination of the 
nondestructive data for these plots showed that a few very large 
trees were heavily influencing the LAI estimate. As an example. 
Plot #21 contained a tree with a DBH of 45 cm - considerably 
larger than the greatest DBH (35.4) in the sample of sacrificed 
trees, thus resulting in an extrapolated (and probably spurious) 
leaf area estimate which contributed about 30% of the entire 
plot's leaf area estimate! 

The effect of estimating the leaf areas of the sacrificed 
trees by sampling and ratioing manifests itself in the v Qk of 
Section 3.2.1. By setting v Qk equal to zero in the evaluation of 
W in Q 2 ^ ^ , one can see what the MSPE of the plot LAI estimates 
would have been had the leaf areas of the sacrificed trees been 
known. On the average, the standard errors (SE) in Table 4, 
would have been reduced by about 7% for the various plots if the 
VQk were zero. The smallest reduction was on Plot#98 (2%), while 
the largest reduction was 13% on Plot 71. 

5.0 CONCLUSIONS AND RECOMMENDATIONS 
In this paper, we have shown how various components of error 
may be identified and quantified in a typical ecological 
estimation problem. While it is true that many of the models 
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employed were approximate, it is felt that they represented 
reality well enough to provide a viable means of assessing the 
relative contributions of the many sampling and estimation errors 
encountered, as long as extrapolation was avoided. The results 
of this analysis can be used to design future surveys of this 
type more efficiently in terms of both sampling strategy and 
determination of sample size. Below, we list some improvements 
that could be made to the COVER experiment that have been 
identified as a direct result of this analysis. 

1. Number of Sacrificed Trees - Given the same sampling 
strategy on the plots, the number of sacrificed trees should be 
at least doubled. Doing so for this study would have reduced the 
number of plots from 20 to 11 for which the dimension analysis 
coefficient estimation error component was dominant. Tripling 
the number of sacrificed trees would leave only five such plots. 

2. Branch Sampling and Regression Within Sacrificed Trees - 
Care should be taken to include all non-representative and 
particularly the largest branches in the sample for which foliage 
weights were measured. The range of the predictor variable (e.g. 
branch depth of crown) for the unsampled branches should be 
well-covered by the values of this variable for the sampled 
branches in the regression. 

3. Area-to-Weight Ratios - Multiple samples of foliage 
should be taken from each crown stratum to provide a means of 
corroborating the additive model (3.6) and for variance 
estimation. At least rough estimates of sampling fractions for 
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the ratios. , 

4. Dimension Analysis - With a larger number of sacrificed 
trees, it would be better to fit separate models for small and 
large trees. For the few extremely large plot trees which are 
bigger than the largest sacrificed trees, extra effort should be 
made to estimate their leaf area by intense non-destructive 
measurements to avoid the sort of extrapolation that occured on 
Plot #21. For example, even a rough assessment such as "twice 
the leaf area of a neighboring smaller tree" would probably be 
more accurate than a gross extrapolation of the fitted dimension 
analysis model. 
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TABLE 1. Sacrificed Tree Data 


Tree No. 

DBH 

H 

N 

NS 

n 

11 

0.90 

2.20 

9 

9 

8 

8 

1.20 

2.80 

18 

14 

13 

10 

1.45 

3.20 

27 

17 

17 

7 

1.80 

3.78 

37 

16 

15 

6 

2.00 

4.60 

34 

16 

15 

9 

2.20 

3.10 

17 

13 

12 

31 

3.40 

5.35 

28 

16 

16 

32 

3.40 

5.70 

25 

19 

18 

30 

3.50 

5.35 

28 

19 

18 

17 

7.30 

9.20 

52 

20 

18 

16 

9.10 

9.40 

35 

19 

17 

15 

10.50 

11.50 

51 

20 

18 

29 

13.00 

16.10 

20 

13 

12 

28 

13.70 

15.90 

20 

13 

12 

13 

15.10 

16.70 

49 

18 

16 

27 

15.40 

17.40 

24 

14 

- 13 

3 

15.80 

15.60 

63 

63 

59 

5 

17.30 

15.50 

67 

67 

62 

1 

19.40 

23.00 

40 

40 

37 

21 

19.50 

19.35 

20 

13 

12 

18 

21.50 

23.10 

66 

15 

12 

24 

22.50 

22.50 

59 

13 

11 

4 

22.60 

18.10 

60 

59 

55 

25 

22.80 

22.40 

33 

12 

11 

2 

23.00 

22.50 

47 

47 

43 

19 

25.10 

23.80 

27 

14 

13 

26 

25.20 

22.50 

49 

10 

9 

12 

27.80 

23.50 

42 

15 

13 

22 

30.20 

23.50 

40 

16 

15 

20 

32.10 

23.80 

64 

13 

12 

14 

32.40 

23.50 

36 

13 

12 

. 23 

35.40 

22.50 

48 

13 

12 
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TABLE 2. Within-tree Regressions: Foliage Weight = b 0 + bjdc + b 2 dc 2 


TREE DBH 

P 1 

b d 

COEFFICIENTS 
b l b 2 

1 

b n 

t 

b 1 

VALUES 

11 

0.90 

2 1 

0.13 

— 

35.53 


0.12 


2 ' 

8 

1.20 

2 1 

0.47 

3.22 

— 


1.86 

1.98 


10 

1.45 

2 I 

0.56 

6.04 

— 


1.06 

2.27 


7 

1.80 

3 1 

0.35 

3.54 

8.78 


0.38 

0 . 50 

0 . E 

6 

2.00 

2 1 

1.26 

— 

19.59 


7.60 


7.3 

9 

2.20‘ 

2 1 

4.06 

— 

23.87 


2.62 


1 C 

31 

3.40 

2 ! 

1.55 

— 

30.48 


3.26 


7 8 

32 

3.40 

2 1 

0.93 

— 

22.83 


2.79 


9 6 

30 

3 . 50 

2 1 

A 

4.64 

11.39 



1.88 

2 8 

17 

7.30 

1 

— — 

37.41 

14.03 


— 

1.84 

0.6 

16 

9.10 

2 1 

— 

28.79 

36.97 



2 . 57 

2 2 

15 

10.50 

2 ! 

— 

33.79 

49.06 




1.71 

2 4 

29 

13.00 

2 1 

— 

49.69 

54.58 




1 .76 

2 1 

28 

13.70 

2 1 

— 

29.31 

91.48 


_ __ 

1.64 

5 1 

13 

15.10 

3 1 

0.50 

47.47 

29.72 


0.40 

4.84 

3.9 

27 

15.40 

3 1 

0.32 

. 59.62 

56.82 


0.11 

2.61 

3 . 9 

3 

15.80 

2 1 

— 

31.64 

34.17 



6.23 

7 3 

5 

17.30 

2 1 

— 

32.18 

22.59 


— 

6.55 

6.2 

1 

19.40 

2 1 

— 

43.97 

28.60 




3.72 

4 . 51 

21 

19.50 

2 1 

— 

36.23 

80.81 


— 

1.15 

4 . 2 : 

18 

21.50 

2 1 

— 

62.01 

33.31 



2.33 

1 8! 

24 

22.50 

3 i 

1.72 

- 12.38 

113.67 


0.54 

0.31 

5. If 

4 

22.60 

3 1 

0.93 

27.61 

34.62 


0.74 

2.75 

6.1" 

25 

22.80 

3 1 

0.71 

54.08 

78.40 


0.36 

2.07 

7 . 3f 

2 

23.00 

2 1 

— 

26.29 

37.82 



2.15 

5.9! 

19 

25.10 

1 | 

— 

140.56 





9.10 


2b 

25.20 

2 1 

— 

50.61 

20.55 



2.14 

1. 5 C 

12 

27.80 

2 1 

— 

28.81 

42.86 




0 .56 

2 9( 

22 

30.20 

2 1 

4.14 

— 

89.18 


0.74 


10.44 

20 

32.10 

3 1 

2.47 

61.71 

46.91 


0.60 

1.15 

2.21 

14 

32.40 

2 1 

— 

16.16 

71.00 



0.74 

7 . 7 C 

23 

35.40 

2 1 


77.40 

82.64 


— 

1.95 

5.81 


36 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Ta£>le 3. Within-Tree Leaf Area Estimation: Error Analysis 


Tree 

DBH 

»’ ' 

A • 

SE 

PROP MS PE 
BRANCH 
SAMPLING 

PROP MSPE 
RATIO 
EST 

11 

0.90 

0.43 

0.05 

.000 

1.000 

8 

1.20 

0.18 

0.03 

.449 

.551 

10 

1.45 

0.37 

0.05 

.460 

.540 

7 

1.80 

0.91 

0.17 

.793 

.207 

6 

2.00 

0.85 

0.10 

.204 

.796 

9 

2.20 

1.12 

0.22 

.737 

.263 

31 

3.40 

3.19 

0.42 

.435 

.565 

32 

3.40 

2.03 

0.15 

.225 

.775 

30 

3.50 

1.41 

0.11 

.435 

.565 

17 

7.30 

10.41 

1.88 

.875 

.125 

16 

9.10 

8.31 

1.15 

.583 

.417 

15 

10.50 

14.32 

1.47 

.412 

.588- 

29 

13.00 

11.01 

1.28 

.772 

.228 

28 

13.70 

10.97 

1.23 

.270 

.731 

13 

15.10 

8.79 

0.82 

.252 

.748 

27 

15.40 

13.94 

1.00 

.448 

.552 

3 

15.80 

19.39 

1.55 

.000 

1.000 

5 

17.30 

21.44 

1.61 

.000 

1.000 

1 

19.40 

31.44 

2.24 

.000 

1.000 

21 

19.50 

17.46 

1.53 

.260 

.740 

18 

21.50 

18.38 

2.24 

.888 

.112 

24 

22.50 

49.93 

5.53 

.487 

.513 

4 

22.60 

28.71 

2.06 

.000 

1.000 

25 

22.80 

41.60 

3.92 

.212 

.788 

2 

23.00 

38.67 

2.49 

.000 

1.000 

19 

25.10 

27.20 

2.85 

.463 

.538 

26 

25.20 

23.71 

4.86 

.818 

.182 

12 

27.80 

72.29 

7.95 

.668 

.333 

22 

30.20 

74.20 

8.35 

.352 

.648 

20 

32.10 

52.49 

8*01 

.795 

.205 

14 

32.40 

102.01 

10.75 

.443 

.557 

23 

35.40 

120.80 

13.29 

.416 

.584 
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TABLE 4 . Plot-Level LAI Estimates and Error Analysis 


: PLOT 

NO. OF ^ 
TREES f 

A 

LAI 

SE 

CV 

SAMP 

VAR 

COFF EST 
VAR 

MODEL 

VAR 

RJ 

99 

319 

0.0889 

1.293 

0.193 

0.1493 

0.0135 

0.0218 

0.0020 

0 . ( 

86 

301 

0.0678 

1.593 

0.315 

0.1977 

0.0626 

0.0331 

0.0032 

l.E 

94 

458 

0.0889 

1.642 

0.227 

0.1382 

0.0140 

0.0352 

0.0024 

o.: 

3 6 

236 

0.3556 

2.047 

0.225 

0.1099 

0.0321 

0.0157 

0.0030 

2 . C 

95 

569 

0.0889 

2.058 

0.439 

0.2133 

0.1344 

0.0553 

0.0031 

2.4 

83 

175 

0.3556 

2.323 

0.211 

0.0908 

0.0204 

0.0206 

0.0034 

0.9 

16 

100 

0.3556 

2.413 

0.223 

0.0924 

0.0269 

0.0190 

0.0038 

1.4 

20 

231 

0.3556 

2.416 

0.429 

0.1776 

0.1507 

0.0304 

0.0031 

4.9 

89 

297 

0.0756 

2.429 

0.598 

0.2462 

0.2750 

0.0769 

0.0058 

3.5 

80 

369 

0.3556 

2.451 

0.240 

0.0979 

0.0337 

0.0202 

0.0037 

1.6 

82 

202 

0.3556 

2.496 

0.255 

0.1022 

0.0281 

0.0336 

0.0033 

0.8 

3 

95 

0.3556 

2.507 

0.244 

0.0973 

0.0362 

0.0193 

0.0041 

1.8 

93 

86 

0.3556 

2.536 

0.246 

0.0970 

0.0349 

0.0215 

0.0039 

1.6; 

87 

194 

0.0656 

2.591 

0.431 

0.1663 

0.0907 

0.0872 

0.0080 

1.0< 

84 

644 

0.0733 

2.634 

0.403 

0.1530 

0.0669 

0.0906 

0.0048 

0.7: 

69 

377 

0.0889 

2.789 

0.406 

0.1456 

0.0575 

0.1013 

0.0059 

0 . 5* 

81 

145 

0.3556 

2.815 

0.212 

0.0753 

0.0135 

0.0271 

0.0043 

0.45 

85 

230 

0.3556 

2.815 

0.294 

0.1044 

0.0137 

0.0699 

0.0031 

0.1S 

71 

498 

0.0889 

2.835 

0.354 

0.1249 

0.0148 

0.1048 

0.0054 

0.14 

75 

16l 

0.3556 

2.911 

0.288 

0.0989 

0.0200 

0.0592 

0.0036 

0.33 

98 

262 

0.3556 

2.924 

0.370 

0.1265 

0.1041 

0.0281 

0.0045 

3.70 

72 

167 

0.3556 

2.933 

0.297 

0.1013 

0.0243 

0.0601 

0.0036 

0.40 

96 

173 

0.3556 

3.033 

0.253 

0.0834 

0.0146 

0.0454 

0.0041 

0.32 

73 

107 

0.3556 

3.034 

0.195 

0.0643 

0.0035 

0.0299 

0.0048 

0.11 

92 

108 

0.3556 

3.071 * 

0.236 

0.0768 

0.0150 

0.0361 

0.0045 

0.41 

21 

144 

0.3556 

3.082 

0.799 

0.2592 

0.5574 

0.0724 

0.0090 

7.69 

74 

142 

0.3556 

3.134 

0.238 

0.0759 

0.0088 

0.0437 

0.0044 

0.20 

90 

186 

0.3556 

3.167 

0.248 

0.0783 

0.0153 

0.0416 

0.0046 

0.36 

97 

16 / 

0.3556 

3.264 

0.305 

0.0934 

0.0247 

0.0642 

0.0042 

0.38 

77 

105 

0.3556 

3.322 

0.263 

0.0792 

0.0297 

0.0339 

0.0055 

0.87 

88 

370 

0.0811 

3.378 

0.533 

0.1578 

0.1279 

0.1482 

0.0083 

0.86: 

79 

283 

0.3556 

3.945 

0.291 

0.0738 

0.0298 

0.0480 

0.0067 

0.621 


•*, ''RATIO" = (SAMP VAR)/(C0EFF EST VAR) 
* 
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Figure 1: Plot Spectral Measurements 
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Figure 2: Tree Dimensions 
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Figure 3: Branch Dimensions 
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Figure 5: Comparison of smoothed to raw ratios: Leaf area to leaf weight 
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Figure 6: Dimension analysis model 
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